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Abstract 

We present closed-form analytic solutions to non-secular Bloch-Redfield master equations for 
quantum dynamics of a V-type system driven by weak coupling to a thermal bath. We focus on 
noise-induced Fano coherences among the excited states induced by incoherent driving of the V- 
system initially in the ground state. For suddenly turned-on incoherent driving, the time evolution 
of the coherences is determined by the damping parameter C, = ^(71 -|- 72 )/Ap, where 7 * are 
the radiative decay rates of the excited levels i = 1 , 2 , and Ap = depends 

on the excited-state level splitting A > 0 and the angle between the transition dipole moments 
in the energy basis. The coherences oscillate as a function of time in the underdamped limit 
(C ;$> 1 ), approach a long-lived quasi-steady state in the overdamped limit (C ^ 1 ), and display 
an intermediate behavior at critical damping = 1). The sudden incoherent turn-on generates 
a mixture of excited eigenstates |ei) and |e 2 ) and their in-phase coherent superposition \f+) = 
-|- y/r 2 \e 2 )), which is remarkably long-lived in the overdamped limit (where ri and r 2 
are the incoherent pumping rates). Formation of this coherent superposition enhances the decay 
rate from the excited states to the ground state. In the strongly asymmetric V-system where the 
coupling strengths between the ground state and the excited states differ significantly, we identify 
additional asymptotic quasistationary coherences, which arise due to slow equilibration of one of 
the excited states. Finally, we demonstrate that noise-induced Fano coherences are maximized with 
respect to populations when ri = r 2 and the transition dipole moments are fully aligned. 
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I. INTRODUCTION 


Weak-field thermal excitation of multilevel quantum systems is the primary step in many 
physical, chemical, and biological phenomena, ranging from photosynthesis [1, 2] to photo¬ 
voltaic energy conversion [3], visual phototransduction [4] and nanoscale heat transfer [5]. 
Several of these phenomena involve absorption of sunlight, followed by conversion of solar 
energy into (electro)chemical energy. Since the primary light reactions steps of photosyn¬ 
thesis are extremely efficient, detailed understanding of the physics of solar energy capture 
and transformation can lead to new insights into the bio-inspired design of highly efficient 
artificial photovoltaic devices (such as solar cells and photodetectors) [6-9]. For this reason, 
this work discusses sunlight excitation. However, our results are relevant to excitation by 
any incoherent source (e.g. incoherent phonons). 

The primary absorption of sunlight in biological systems occurs in photosynthetic light¬ 
harvesting complexes (LHCs) which also funnel the excitation energy to the reaction center, 
where it is used to drive the subsequent "dark" reactions of photosynthesis [1, 2, 10]. The 
initial stages of photosynthetic light harvesting can be probed by initiating and measuring 
the dynamics of molecular excitations in real time via advanced spectroscopic techniques 
such as the two-dimensional photon echo spectroscopy (2DPE) [2, 16]. Recent 2DPE ex¬ 
periments revealed long-lasting, wavelike dynamics of energy transfer between the different 
chromophores in the Fenna-Matthews-Olson (FMO) complex of green sulphur bacteria [12] 
and in the PC645 complex of photosynthetic algae [13]. These room-temperature 2DPE 
experiments yielded coherence lifetimes in excess of 500 fs, much longer than expected for 
the electronic coherences in typical condensed-phase environments at room temperature 
[2]. While these surprising discoveries stimulated a wave of theoretical and experimental 
research into quantum coherent effects in LHCs [2, 14, 15], the origin of the observed coher¬ 
ences (electronic vs. vibrational) and their relevance to the naturally occurring biological 
processes remains a subject of vigorous debate [16]. 

The concern regarding the role of quantum coherences in solar light-harvesting is relevant 
because 2DPE experiments use ultrafast transform-limited light pulses to excite coherent, 
time-evolving superpositions of molecular eigenstates (wavepackets). In contrast, the Sun 
emits blackbody radiation at T = 5800 K, which is a statistical mixture of number states [17], 
characterized by a diagonal density matrix in the number-state representation [18]. Because 
of the absence of a phase relation between the different components of the mixture, sunlight 
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is expected to populate molecular eigenstates, not producing any coherences between them. 
Using these arguments, which are in keeping with the standard Einstein theory of light- 
matter interaction [23], several groups showed that there is no dynamical evolution of the 
excited states prepared by incoherent light [19, 20]. More recent results [21] do suggest, 
however, possible time dependence at shorter times in complex molecular systems. 

The most general theory of weak-field incoherent excitation of multilevel quantum sys¬ 
tems is based on the Bloch-Redfield (BR) quantum master equations, which describe the 
time evolution of the reduced density operator of a few-level system interacting with a 
continuum of quantized radiation field modes [22]. The Pauli rate equations underlying 
the Einstein theory [23] follow from these equations upon neglecting the non-secular terms, 
which describe the coupling between the populations and coherences. The non-secular terms 
are responsible for the phenomenon of Fano interference between the different incoherent ex¬ 
citation pathways leading to the same final states [24, 25, 27, 28]. Similar coherences (albeit 
induced by laser excitation) give rise to a host of remarkable phenomena in quantum optics 
such as electromagnetically induced transparency [29, 30], slow light [29, 31], and lasing 
without inversion [30, 32]. The existence of Fano coherences in few-level quantum optical 
systems driven by incoherent light [24, 25, 27] naturally raises the question of their rele¬ 
vance in the context of natural light-harvesting and artificial photovoltaics, where they have 
been proposed to enhance the efficiency of quantum heat engines [8, 9]. However, despite 
early recognition of these noise-induced coherences [24], their properties remain incompletely 
understood. 

Previous theoretical work has focused on steady-state Fano coherences in model V- (or 
A)-type systems with degenerate upper (or lower) levels [27]. In particular, Kozlov et al. 
demonstrated that incoherent excitation can lead to non-vanishing steady-state coherence 
between the degenerate upper levels [25]. Agarwal and Menon [36] showed that the noise- 
induced Fano coherences between non-degenerate eigenstates vanish in the long-time limit, 
leading to a canonical steady state. They also identified the signatures of the coherences in 
steady-state ffuorescence spectra. Hegerfeldt and Plenio demonstrated that Fano coherences 
can give rise to quantum beats and to an extended dark period in time-resolved fluorescence 
emitted by V- and A-type systems excited by incoherent light [37], while Altenmiiller clarified 
the origin of these quantum beats as originating from a factorized system-reservoir initial 
state [38]. Similar results have been obtained for the A-type system [27, 35, 39], where the 
coherence can be generated by spontaneous decay to closely spaced ground levels. 
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We have recently provided analytical solutions to the BR equations for an arbitrary 
excited-state spacing A and identified two important regimes of incoherent excitation for 
the specific case of equal pumping rates [40]. In the large-spacing regime A/y ^ 1, the 
coherences between the excited eigenstates display an oscillating behavior, decaying on the 
timescale = I/ 7 , where 7 is the spontaneous decay rate of the excited states. In the 
opposite regime of closely spaced excited levels A /7 1 , the coherences decay on the 

timescale ta = 7 /( 2 A^), thus suggesting that near-degenerate energy levels can maintain 
coherence on a very long timescale [40]. However, this preliminary work did not consider 
several factors that affect the coherence dynamics: ( 1 ) the relative magnitude of incoher¬ 
ent pumping rates, and ( 2 ) the magnitude of the population-to-coherence coupling (or the 
transition dipole moment alignment factor p = ■ ^ige 2 /(lh 3 eiIIhge 2 !)• Here are the 

transition dipole moments between the ground and excited state [25, 27, 40]. While it 
is known that the Fano coherence gets smaller with decreasing p and eventually vanishes at 
p = 0, the exact dependence on p remains unexplored, as is the case of unequal pumping 
rates. Investigating noise-induced coherent dynamics away from the p = 1 limit explored 
previously [ 8 , 9, 25, 40] is particularly important for molecular systems, where transition 
dipole moments will typically have a broad range of magnitudes and orientations. 

Below, we present analytic solutions to the BR equations of motion describing incoherent 
excitation of the V-system for arbitrary values of pumping rates and p. Extending our 
previous work [40], we find that in the weak-pump limit the solutions can be classified in 
terms of two parameters defined in analogy to the damped harmonic oscillator: ( 1 ) (, which 
characterizes the dynamics of the excited state coherences in direct analogy to the damping 
coefficient of a classical oscillator, and ( 2 ) 7 , which quantifies the anharmonicity of the 
coherent dynamics. We show that the relevant parameter, which governs the dynamical 
evolution of noise-induced coherences in the asymmetric V-system is the ratio = Ap /7 
of “renormalized” excited-state splitting Ap = -y/A^ -(- (1 — p^) 7 i 72 to the average decay 
width 7 . We identify two key dynamical regimes: underdamped ((/ <^ 1), where the excited- 
state coherences oscillate multiple times before decaying to zero on the timescale = 
1 / 7 , and overdamped {( <C 1 ), where the coherences approach a quasisteady state after a 
short time interval (r^) following the initial turn-on of incoherent driving, and then stay 
constant for Tiong = l/^p- In this regime, coherences can be made arbitrarily long-lived 
by reducing Ap. A critical damping regime occurs at (/ = 1 , where the coherences display 
an intermediate behavior between the overdamped and underdamped regimes. We present 
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analytic expressions for the noise-induced coherences in each of these limits, and find that the 
coherences are maximized when the pumping rates are equal to each other. Taken together, 
these results completely characterize the dynamical regimes of the V-system weakly driven 
by suddenly turned on incoherent light. (Issues associated with slow turn-on are under 
consideration [33].) 

We note that this work provides new physical insights into the dynamics of the V-system 
driven by incoherent radiation, which go well beyond those reported previously (including 
Ref. [40]). In addition to providing novel analytic expressions for density matrix dynam¬ 
ics valid for arbitrary pumping rates, transition dipole orientations, and level splittings 
[Eqs. (10)-(16)], this Article 

1. Establishes the existence of three important dynamical regimes: underdamped, criti¬ 
cal, and overdamped, classifiable by the value of the damping parameter ( in Eq. (8b). 

2. Explores the dynamics of the critically damped V-system (Figs. 7, 8 and Sec. IIIC [34]). 

3. Establishes the presence of asymptotic quasi-stationary coherence in the asymmetric 
V-system (Sec. IIIA and Figs. 3E-F). 

4. Provides a physical interpretation of the quantum state created by incoherent pumping 
as a mixture of excited eigenstates and their coherent superposition [Eq. (3)]. 

5. Demonstrates that the noise-induced coherences are maximized when the pumping 
rates are equal and the transition dipoles are fully aligned. 

This paper is organized as follows. Section II presents the BR equations of motion, 
outlines some general properties of the system dynamics, and gives our general weak-field 
solution for the V-type system driven by an incoherent radiation field. Closed-form analytical 
solutions for the time dynamics of populations and coherences in the limiting case of highly 
separated excited states ((^ 1, valid for small molecules) are presented in Sec. IIIA. 

Section IIIB focuses on the case of 3> 1, which holds for large molecules with closely 
spaced vibronic energy levels. Section IIIC considers a special intermediate regime where 
the spacing of the excited states is comparable to their radiative decay widths {( = 1). The 
results for varying alignment parameter p and the pumping rates are presented in these 
subsections. Section IV concludes with a brief summary of main results and outlines several 
directions for future work. 
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Fig. 1: Schematic representation of a V-type System. A is the excited state splitting, 7 * is the 
radiative line-width, and rj is the incoherent pumping rate of excited state |ej). 

II. BLOCH-REDFIELD MASTER EQUATIONS 

Consider a V-system interacting with a suddenly turned-on incoherent radiation held, 
under the Born-Markov approximation. The dynamics of such a system is governed by the 
following system of BR master equations [40, 41] 

Pe,e^ = -(u + + UPgg - p{y/n^ + ^/^)pe,e2 (1^) 

1 p p 

Peie2 = - 2 (^1 +^2 +7l +72)Peie2 - A + -y/¥^{2pgg-pe^ei -Pe2e2 ) “ ^ +^6262 ) 

(lb) 

where here and below atomic {h = 1) units are used. In Eq. (1), absorption and stimulated 
emission processes are parametrized by the incoherent pumping rates of the \g) -H- |ej) transi¬ 
tions Ti = BiW{u!ge^), given by the product of the Einstein R-coefhcients Bi = 7i\pge^\^/{3eo) 
and the intensity of the incident blackbody radiation at the corresponding transition fre¬ 
quencies W{(JJge^)■ Spontaneous emission processes are governed by the radiative decay 
widths of the excited states, 7* = a;ggJ/igeiP/( 37 reoc^), A = a;eie2 gives the excited state 
energy splitting, and p = /j.gei ■ P-ge2li\pgei\\Pge2\) quantihes the alignment of the \g) -H- \ei) 
transition dipole moments, Rgei- The notation denotes the real part of Peie2 while 

Peie2 denotes the imaginary part. Here we neglect the environment-induced dephasing and 
relaxation processes, assuming that the rates of excited state relaxation and dephasing are 
small compared to those of the radiative processes (absorption, decay and stimulated emis¬ 
sion). Eq. (1) makes the usual assumptions of weak system-bath coupling and short (delta 
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function) time-correlation in the radiation field [22, 32], 

The one-photon coherences between the ground and excited states oscillate extremely 
quickly and can be decoupled from the dynamics of the rest of the system (the partial 
secular approximation) [40, 41], The quantum master equations [Eq. la] for the populations 
contain secular rate-law terms which depend only on the diagonal density matrix elements, 
PeieiJ coherence-dependent non-secular terms (proportional to Peie 2 ) that lead to the 
breaking of detailed balance responsible for the proposed efficiency enhancement in quantum 
heat engines [ 8 , 9]. The rate-law terms contain the effects of the independent \g) -H- \ei) 
transitions while the coherence-dependent terms show the effect of interference between the 
two pathways on the emission processes. The dynamics of the V-system can be understood 
in all regimes in terms of the interplay of interference and rate-law transitions. 

The quantum master equations can be recast in the Liouville representation in terms of 
the state vector x = [Peiei, Pezez, pfjea’ where are the real and imaginary 

parts of the excited state coherence, respectively. Substituting the normalization condition 
[pgg = 1 — Peiei — Pe 2 e 2 ) Eq. (1) allows US to Write the master equations in the following 
inhomogeneous form; 


A 


-(n + 7i) 
-r2 

(3^^+ 


—X = Ax -F d 
at 


-ri 

-(^^2 + 72 ) 
( 3 ^/f 7 r^ + V7172) 


-p + y/frn) 

-P + V7172) 

-(f+ 7 ) 


( 2 a) 

0 ^ 

0 

A 


v 


0 


0 

/ A \ 


-A -(?^ + 7 ) / 

( 2 b) 


d = 


r2 

Py/rir2 


( 2 c) 


V 0 

where f = ^(ri-|-r 2 ) and 7 = \{'^i+'^ 2 ) are the mean pumping rate and radiative decay width, 
respectively. In Eq. (2), A is a coefficient matrix, analogous to the Liouville superoperator, 
and d is a driving vector containing the contributions of radiative excitation from the ground 
state Ip). The driving vector, d (Eq. 2c), indicates that the incoherent radiation drives the 
V-system into the following mixed state on the excited state manifold: 


Pd oc (1 -p)(ri|ei)(ei| + r 2 |ei)(e 2 |) + p|0+)((/>+| (3) 
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where |0+) = :^(\/n|ei) + v^|e2)) is an in-phase coherent superposition of the energy 
eigenstates. Equivalently, only the real part of the coherences, directly excited from 

the ground state [see Eq. (lb)]. Equation (3) illustrates the important role of the alignment 
factor p in parametrizing the preparation of the coherent superposition |0+). In contrast, 
the Pauli rate-law equations predict excitation into the incoherent mixture of excited states: 

Pincoh oc ri|ei)(ei| -1- r 2 |e 2 )(e 2 | (4) 

that equilibrates to the mixture: 


Peq OC |ei)(ei| -f |e2)(e2| 


(5) 


Physically, |0+) is generated by a simultaneous excitation from the ground state to the two 
excited states |ei) and |e2). In the limit of orthogonal transition dipole moments {p = 0), 
Eq. (3) gives an incoherent mixture of the energy eigenstates [Eq. (4)] while in the opposite 
limit of parallel dipole moments {p = 1), the pure in-phase coherent superposition, |0_|_)((/)+|, 
is produced. 

Using the standard variation of parameters procedure [42], the solutions to the master 
equations (Eq. (2c)) with initial conditions corresponding to excitation from the ground 
state (x(t = 0) = 0) is given by 

x(t) = [ ^ V / ds{v, ■ d)e^di-^)v^ (6) 

Jo Jo 

where Aj is the eigenvalue of A and Vj the corresponding eigenvector. Eq. 6 relates the 
eigenvalues, {Aj}, of A to the timescales of the system’s evolution r* = Re(Aj)“^ and the 
frequencies of its oscillations oji = Im(Ai). 

In general, the spectrum of the coefficient matrix. A, (Eq. 2b) is complicated. However, 
we note that the intensity of the incident blackbody radiation encountered in nature (rep¬ 
resented here by the effective thermal occupation number ft = —), is typically very weak, 

'yi ' 

ft 1 (the weak pumping limit). This allows us to neglect the terms of order r* in Eq. (2b) 
and greatly simplifies the calculation of the spectrum, {Aj} of A, to give: 


A. 



AJ)±^(V-A2)^ + A^(7,-7.)2 
2 


(7) 


where Ap = + (1 — p‘^)jij 2 and 7 = ^(71 + 72)- The algebraic details of the derivation 

of Eq. (7) are given in Eq. (Al) of the Appendix. 







The evolution of the coherences, Peie2) can be understood in analogy to the behavior of a 
damped harmonic oscillator. To illustrate this analogy, consider a system with symmetrically 
broadened excited states (71 = 72). In this case, Eq. (7) gives one degenerate negative root 


(Ai = —7) and a pair of (possibly complex) roots (A2,3 = —7 ± y (7^ — A^)) analogous to a 
damped oscillator with natural frequency Ap and damping constant 7 [43]. 

Motivated by the similarities to the damped oscillator, we rewrite Eq. (7) as 



(8a) 



(8b) 


A|7 i - 72I 


(8c) 



where the parameter ( plays the same role as the damping coefficient of a classical damped 
oscillator. On the other hand, Eq. (8c) defines a second parameter, 7, which quantifies the 
deviation from the damped oscillator spectrum, (with rj = 0 corresponding to the symmet¬ 
rically pumped V-system, which is equivalent to the damped harmonic oscillator). As such, 
Eq. (Eq. (8c)) relates the asymmetry of the excited states (in terms of their splitting and 
radiative broadening) to their anharmonicity, rj. 

To clarify the effects of rj on the spectrnm Eq. (Eq. (8a)) and hence its deviation from 
that of a harmoni c oscillator consider the positive and negative branches of the rj dependent 
term ('\/^^^^^) separately. In the underdamped regime, the negative branch in Eq. (8a) 
splits the degenerate eigenvalne to give two new relaxation rates, one slower than the original 
rate and the other faster, while the positive branch increases the frequency of the system’s 
oscillations. On the other hand, the negative branch (1 — \Jl + if') leads to asymmetry 
induced oscillations in the overdamped regime S> 1). The positive branch modifies the 
relaxation rates in this regime, increasing the faster rate while slowing down the slower rate. 

To separate the effects of the asymmetry of 71 and 72 from their magnitude, 27 we 
introduce the parameter f3 such that 


7i = 27sin2(/3) 


(9a) 


72 = 27cos^(/3) 
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(9b) 

(9c) 













For fixed p, f5 provides a convenient independent measure of the asymmetry of the transitions, 
with decreasing f3 G [0, indicating a larger difference in decay rates I 71 — 72 |- 

In previously reported analytical results on the V-system, derived for symmetrically 
broadened states (71 = 7 = 72) and collinear dipole moments [p = 1 ), the parameter 
A/7 was introduced [40], analogous to = Ap/7. However, to highlight the similarities to 
the damped oscillator we opt for the parameter (. Figure 2 sketches the various dynamical 
regimes of the V-system excited by incoherent light. The vertical line corresponds to (/ = 1 
and separates the underdamped and underdamped regimes. The region below the horizontal 
line /3 = 7 r /4 corresponds to the asymmetric V-system, with the degree of anharmonicity 
increasing as [3 decreases. As discussed below, the coherent dynamics in the strongly an- 
harmonic region (corresponding to ^ 0 ) displays asymmetry-induced quasistationary 
coherences that are not present in the symmetric V-system. 


III. ANALYTICAL SOLUTIONS 
A. Under damped Regime C ^ 1 

Consider hrst a weakly damped system with C ^ 1 (or equivalently ^ S> 1). Physically, 
this corresponds to a system with a large excited state splitting. A, such as a small molecule 
with few vibronic degrees of freedom. In this limit, the spectrum, {A*}, and corresponding 
normal modes, {vj}, of A can be calculated through a binomial expansion of Eq. ( 8 a) and 
substituted into Eq. ( 6 ) to give the density matrix elements. (The details of this calculation 
are provided in section Eq. (A 2) of the Appendix.) This gives 

PeueAt) = ( 10 a) 

Peue2A) = e~^*sin(Apt) ( 10 b) 

L\p 

^ - 1 ) - £ 11 ^) ( 100 ) 

where he is the effective thermal occupation {e.g. he = a[e'^BT — 1 ]“^ = an. Here a is 
some scalar and h is the thermal occupation of electromagnetic held modes at the ground- 
to-excited manifold transition frequency, cjq (see Fig. 1). We have assumed, without loss of 
generality, that 71 <72. We use the same occupation number for the \g) -H- |ei) and \g) -H- 
1 62 ) transitions. This is justihed since A <C Wq under the Wigner-Weisskopf approximation 
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Fig. 2: A "phase diagram" depicting the regions in parameter space corresponding to the regimes 
of the V-system. The underdamped regime C < 1 shows damped oscillatory coherences, while the 
overdamped and critically damped regimes C ^ 1 show no oscillations. As in Eq. (8b), C = 

/3 is defined by 71 = 27 cos^(/ 3 ), 71 = 27 sin^(/?) (Eq. 9). Note that the difference in line widths 
(I 71 — 72 I) increases as j3 G [0, decreases. 

[32, 40], leading to = n^'^i or ri /71 = r 2/'^2 = ffe, which guarantees that a coherence-free 
(canonical-like) state of the V-system is reached in the limit f —)■ 00 for non-degenerate 
excited states [36]. 

The evolution of the density matrix elements is illustrated in Fig. 3. In this underdamped 
limit, Eq. (10a) and Fig. 3A-C show the evolution of the excited state populations, Peiei, 
to a steady state value, at a rate 7 *, as predicted by the Pauli rate equations. The lack 
of p dependence in Eq. (10a) and the adherence to the dynamics predicted by the secular 
rate-law equations indicates that the populations of a V-type system in the underdamped 
regime is unaffected by the coherences. In other words, the time dynamics of the diagonal 
elements of the density matrix is described equally well by the Pauli-type rate-law equations 
[i.e. Eq. (1) with p = 0]. 
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By contrast, Eq. (10b) and Figs. 3D-F show damped oscillatory coherences, Peie 2 - Equa¬ 
tions ( 10 b) and ( 10 c) can be understood in terms of the coherent | 0 _|_) = -^(ri|ei) +r 2 \e 2 )) 
component of the mixture, pd, prepared by the incoherent radiation, Eq. (3). The |0+) com¬ 
ponent of pd is initially in an in-phase (positive) superposition of the excited states. The 
unitary evolution induced by the (isolated) V-type system Hamiltonian causes the compo¬ 
nents in |ei) and |e 2 ) to accrue a relative phase with frequency A, producing the oscillations 
seen in Eqs. (10b) and (10c) and in Figs. 3D-F. However, the interaction with the radia¬ 
tion held leads to a decay of the coherent wavepacket through spontaneous and stimulated 
emission, leading to exponential decay at a rate 7 , as seen in Eqs. (10b) and (10c) and 
Figs. 3.D-F. Eventually, the interaction with the radiation held destroys the coherent |(/)+) 
component of the mixture generated by the blackbody radiation pd Eq. (3) to give the 
incoherent mixture predicted by the rate-law equations, peq [Eq. (5) or (1) under the secular 
approximation (p = 0 )]. 

A unique feature of the asymmetric V-system (/3 7 ^ 7 r/ 4 ) is the presence of asymp¬ 
totic quasistationary imaginary coherences (Figs. 3E-F). These are contained in the non- 
oscillatory exponentials in Eq. (10c) and are a manifestation of the slow equilibration of one 
of the states. When one state, say |e 2 ), equilibrates much more quickly than the other (i.e. 
72 ^ 7i) then the large population of |e 2 ), Pe 2 e 2 , suppresses the real part of the coherences, 
through the non-secular term coupling to Pe 2 e 2 iii Eq. (lb). Physically, this 

corresponds to constructive interference in the decay processes of in phase superpositions 
(e.g. 10+)), which are characterized by pf^g^ > 0. Consequently, the real coherences decay 
more quickly than the imaginary coherences which survive for longer times due to the slow 
equilibration of |ei). 

Since the coherences are only present in the |0+) components of pd, the p dependence of 
the amplitude of coherences in Eqs. (10b) and (10c) reflects the amount of |0+) in the pd 
mixture, Eq. (3). Further, the magnitude of coherences in a pure |0+) = +^ 2 |e 2 )) 

state is proportional to This produces the scaling in the amplitude of the 

coherences, peie 2 iii Eqs. (10b) and (10c) and hence the decreased amplitude of coherences 
as the difference between 71 = ^ and 72 = ^ increases in Figs. 3.D-F. 

The coherence ratio C{t) = ^ ^ — provides a convenient method of quantifying 

the magnitude of the coherences relative to the populations of the corresponding states 
[40, 44, 45]. Neglecting the weak asymptotic imaginary coherences that appear in highly 
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asymmetric systems, the coherence ratio for an under damped V-system is given by 


at )« { 


3-7* 


2 / 


1 A II -t I sin" 1 —I (11) 

V Ap J \l- ) \ 2 ' ^ ’ 

and is shown in Fig. 4. The C-ratios exhibit damped oscillations, decreasing at short time 
as l/(Apt), extending our earlier result [44, 45, 45]. 

It follows from Eq. (11) that C{t) is maximized when the decay widths, y*, (and hence 
the pumping rates r^) are equal. Furthermore, it can be seen that the amplitude of the 
coherence ratio scales linearly with the alignment parameter, p, and inversely with excited 
state splitting, A contained in Ap. This scaling presented in the prefactor of Eq. (11), 
can also be seen in Fig. 4. Furthermore, the effects of the asymmetry in y*, that are not 
accounted for in Eq. (11), manifest in Fig. 4C in the small residual value of C{t) between 
Ti = yf^ and T 2 = ys'^ 


B. Overdamped Regime C 1 


Consider now the case of a system that is strongly damped (^ ^ 1 or equivalently 
^ -C 1). Using a similar Taylor expansion of Eq. ( 8 a) as in the underdamped case, the 
evolution of the density matrix is obtained after some algebra, detailed in section Eq. (Al4c) 
of the Appendix, as 


d6i56j(^) ^ — 

2 y 




u(l 

P^/rir2 

27 


—2'yt 


) +c(l 


e 27 *) 


rw* - 


( 12 a) 


(12b) 


where i,j = 1,2 and i ^ j- The imaginary coherence, pfi e 2 ~ ^Peie 2 ) is suppressed by a 
factor -C 1. These density matrix elements are plotted in Fig. 5. 

Equation (12b) and Fig. 5 show that, in the overdamped regime, the V-system first evolves 
to the mixture, pd [Eq. (3)], excited by the incoherent radiation over a time Tshort = (2y)~^. 
Then pd decays to the incoherent mixture, peg [Eq. (5)], generated by the rate-law model. 
However, Eq. (12b) and Fig. 5 show that the coherences (and hence the coherent 10+) 
component of pd) are remarkably long-lived, surviving for 


'^long 


A2 

p 


(13) 


In particular, Eq. (13) indicates that Tiong —)■ cx 3 as Ap 0 and y 7 ^ 0. In the Ap = 0 case, 
Pd is dynamically stable and survives in the long time steady state [25]. Noting that Ap —)■ 0 
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Fig. 3: Evolution of populations and coherences of an underdamped V-system (( <C 1), typical of 
small molecule, evaluated with aligned transition dipole moments (p = 1). Panels A-C show the 
evolution of the populations Peiei (solid) and Pe 2 e 2 (dashed) where state |e 2 ) has the highest energy. 
Subplots D-F show the evolution of the real (p^ solid line) and imaginary g^ dashed line) 
coherences. Here 7 = 1.0 and A = 12.0, with varying pumping rate asymmetry 71 = 7 sin^(/ 3 ) and 
72 = 7cos2(/3). 

indicates that A = 0 and p = 1, this is eqnivalent to the steady state coherences observed 
in previous investigations of the V-system [25, 40], which give rise to coherent population 
trapping and multiple (initial-state-dependent) steady states [25, 26]. 

In a degenerate V-system (i.e., A = 0), the coherent superposition |0_|_) accrues no 
relative phase, corresponding to the vanishing of the coupling between the real and imaginary 
coherences in Eq. (lb). As A —)■ 0, the evolution of the density matrix is determined by the 
interaction between the components of p^. In particular, the eigenstate components, |ei) 
and |e 2 ), suppress the coherent | 0 _|_) component through the population-coherence coupling 
term in Eq. (lb). If p 7 ^ 1 then Eq. (3) indicates that |ei) and |e 2 ) will be present in pd, 
increasing the decoherence rate. However, if p = 1 and A = 0 then the coherent pure state 
Pd = 10+) (0+1 does not decay, producing the steady state coherences observed in Eq. (12b) 
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Fig. 4: Evolution of the coherence ratio C{t) = -—— of an underdamped V-system <C 1). 
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Here 7 = 1.0 and A = 12.0, with varying pumping rate asymmetry 71 = 7 sin^(/ 3 ) and 72 = 
7 C 0 S^(/ 3 ) and dipole alignment, p = 1.0 (solid) and p = 0.7 (dashed). 


|26|, 


The coherence ratio, Cit) = —[AifsJ— Qf V-system in the overdamped regime is given 

Pejej+Pe2e2 


by 


C{t) = 


2f 


e 27 * _ g 274 


a2 

2 — 

A plot of C{t) is shown in Fig. 6 . Both Equation (14) and Fig. 6 illustrate a notably long- 
lived substantial C{t) followed by the collapse of to the incoherent mixture peq in time 
Tiong- Furthermore, note that Eqs. (11) and (14) both show the same scaling dependence on 
Py/^if' 2 - That is, as in the underdamped regime, the p dependence reflects the dependence 
on the amount of | 0 +) present in p^ while gives the coherence of the | 0 +)( 0 +| pure 

state. 


(14) 


C. Critical Regime C = 1? = 0 

Finally, to connect the overdamped and underdamped regions, consider a special case 
between them. In particular, when (^ = ^ = 1 (or equivalently ^) and p = 
the spectrum, {Aj}, of A (Eq. ( 8 a)) becomes completely degenerate with one eigenvalue 7 . 
These constraints are equivalent to 71 = 7 = 72 and A = pp. After some algebra, detailed 
in section Eq. (Al9b) of the appendix, the matrix A can be put into Jordan Canonical Form 
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Fig. 5 : Evolution of populations and coherences of an overdamped V-system ((^ S> 1 ) evaluated 
with aligned transition dipole moment operators {p = 1).Panels A-C show the evolution of the 
populations Peiei (solid) and Pe2e2 (dashed) where state |e2) has the highest energy. Panels D- 
F show the evolution of the real (p^ £2’ solid) and imaginary {pi^ dashed) coherences. Note 
that the imaginary coherences are heavily suppressed and nearly vanish throughout the system’s 
evolution. Here 7 = 0 . 5 , and A = 0 . 0012 , with varying pumping-rate asymmetry 71 = 7sin^(/3) 
and 72 = 7 cos^(/ 3 ) 


and substituted into Eq. (6) to give the density matrix elements. 


Pe,,eAt) -npi e ^ t^e ^‘) 

(15a) 

pfi,e2(^) = nePite~^^ 

(15b) 

PLe2it) = 

(15c) 


The long-time {t —)■ 00 ) limit of Eq. (15c) yields the incoherent mixture of excited states 
predicted by the Pauli rate equations, peg (Eq. (5)). However, in contrast to the under¬ 
damped regime, the coherences alter the approach of the system to its equilibrated state. In 
Eq. (15c), this can be seen in the terms proportional to t and The altered dynamics of the 
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Fig. 6: Evolution of the coherence ratio C{t) = - —— of an overdamped V-system (C !)■ 
Here 7 = 1.0, and A = 0.0012, with varying pumping-rate asymmetry 71 = 7sin^(/3) and 72 = 
7C0S^(/3) and dipole alignment, p = 1.0 (solid) and p = 0.7 (dashed). 

density matrix elements are also clearly seen in Fig. 7. In particular, the slowly equilibrating 
populations, Peiei, show a marked "bend" in their approach to the steady state, resulting 
in a transient suppression of population compared to their underdamped or coherence free 
evolution. More generally, as the damping of the system becomes stronger, the effects of 
the coherences become more prominent and long-lasting. In the critically damped regime, 
no quasistationary behavior is observed. 

In this case, the coherence ratio, C(t) = — jg 71 = 72 

P 7 / e~'^^t^yT^r~A?Wj4: 

Y I 1 - e-^‘ +A 2 t 72 e-^‘ 

The ratio C{t) in the critical region is shown in Fig. 8. Equation (16) shows the same p 
scaling as in the other regimes. However, the dependence on 7 ^ is hidden due to the 71 = 72 
assumption. Figure 8 shows that the same suppression of C{t) with increasing asymmetry 
as in the underdamped (Fig. 4) and overdamped(Fig. 6 ) regions. 




IV. CONCLUSION 

We have presented analytical solutions for the dynamics of a V-system interacting with 
an incoherent radiation bath in three limiting cases. In all three of these regimes the dy¬ 
namics can be understood through the interplay of interference effects and radiative decay 


17 












xio~^A. 0 = 7r/4 


q-B. /3 = 7r/6 


, xio€. 0 = 7^/12 



X10: 


2 

(N 

“ 0 
-1 
-2 


a: 

Ci. 


0 


10 

t (ns) 


20 


10 

t fnsl 


20 



2 


2 

"1 

^-- 


0 

^-1 

-2 


^ 0 

-2 



10 

t (ns) 


20 


Fig. 7 : Evolution of populations and coherences of a critically damped V-system S> 1 ) evaluated 
with aligned transition dipole moment operators {p = 1 ). Panels A-C show the evolution of the 
populations Peiei (solid) and Pe2e2 (dashed) where state |e2) has the highest energy. Panels D-F 
show the evolution of the real (p^ solid) and imaginary (pg^ g^, dashed) coherences. Here 7 = 1 . 0 , 
and A = 1 . 0 , with varying pumping-rate asymmetry 71 = 7sin^(/3) and 72 = 7COS^(/3). 

processes. Our results show a transient suppression of population in the presence of co¬ 
herences relative to rate law dynamics in the overdamped and critically damped regimes. 
These hndings indicate that the in-phase superposition, 10+) prepared in the excited state 
manifold by incoherent light excitation shows constructive interference of radiative decay 
processes, leading to the reduced population of the excited state manifold. Coherent effects 
are transient in all regimes (with the exception of the Ap = 0 case) and the system eventually 
collapses into the incoherent mixtures predicted by Pauli rate-law considerations. However, 
in the overdamped region (0 = 7 /Ap S> 1 ) the coherent superposition | 0 +) is remarkably 
long-lived and survives for Tiong = 27 /(Ap. 

A striking feature of the Bloch-Redheld Master Equations [Eq. (lb)] is that incoherent 
light excitation produces mixtures, pd [Eq. (3)], in the excited state manifold that contain 
in-phase coherent superpositions of the excited states, |0+) = ;^(ri|ei) -1- ^ 2162 )). In- 
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Fig. 8: Evolution of the coherence ratio C{t) = —— q£ critically damped V-system <C 
1 ). Here 7 = 1.0 and A = 1 . 0 , but with varying pumping rate asymmetry 71 = 7sin^(/3) and 
72 = 7 C0S^(/3) and dipole alignment, p = 1.0 (solid) and p = 0.7 (dashed). 

phase superpositions, such as \4>+), lead to constructive interference in the emission process, 
thereby suppressing the population of the excited state manifold as seen in the overdamped 
and critically damped regimes [Eqs. (12a) and (15a) and Figs. 5 and 7]. This observation, 
however, does not contradict the previously reported enhancement of photocell currents due 
to coherences between the excited states [8, 9] as the photocell systems contain an additional 
conduction-band reservoir state |a) to which the two excited states decay, and to which the 
coherent enhancement of decay also applies. The observed enhancement of photocurrent 
in these systems therefore suggests that, rather than enhancing absorption or suppressing 
the radiative decay to the ground state, the coherent enhancement of the decay rate to the 
conduction-band reservoir state la) outweighs the increased decay rate to the ground state. 

Our findings motivate further consideration of the incoherent-light excitation of multilevel 
systems. First, the variation of parameters method used here allows for the consideration 
of time-dependent radiation helds. The investigation of adiabatic turn-on of the radiation 
held [21, 33] is particularly important in the context of natural (sunlight) excitation of 
biomolecules and LHC’s, since the natural turn-on timescales are far longer than typically 
exploited in femtosecond experiments. Second, the Bloch-Redheld method employed here 
predicts rich dynamics in systems with more complicated ground state manifolds, where 
the non-secular terms produce additional phenomena due to interference effects between 
different ground states. Work on these extensions is in progress. 
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Appendix A: Derivation of V System Dynamics 
1. Coefficient Matrix Spectrum 

Consider first the coefficient matrix, A, in the (n 1) weak pumping limit. Equation (2b) 
can be rewritten as a simpler matrix perturbed by a term of order h <C 1. 


/ 

-7i 

0 


0 \ 


^ -71 

-71 

-V^/l\l2 

0 

\ 


0 

-72 

-P\/7i72 

0 

+n 

-72 

-72 

-v^iTn 

0 




-i\/^ 

-7 

A 

-f 

-f ^/^ 

-7 

0 


V 

0 

0 

-A 

- 1 ) 


1 0 

0 

0 

-7 

/ 


(Al) 


where 7 = 1(71 + 72) and f = \{ri + r2) are the arithmetic mean decay widths and pumping 
rate respectively. We treat the perturbative expansion of A given by Eq. (Al) to zeroth 
order in ft +; 1. The results obtained to zeroth order in n are in extremely close agreement 
with the numerically exact solutions and clearly illustrate the physics of the system. 

Consider now the spectrum, {Aj}, of to obtain Eq. (7) in the main text. After some 
elementary manipulations of det(A(°^ — AJ), the characteristic polynomial of takes the 
form 

det(A(°) - AJ) = (71 + A)(72 + A )[(7 + A)^ + A^] -^^ 7172(7 + A)^ (A2) 

Equation (A2) can be rewritten in the following biquadratic form; 

det(A^°^ - XI) = + [ 7 ^ + (1 - P^)lil2 + A^jx + [ 7^(1 - p^)7i72 + 7 i 72 A^] (A3) 

where x = A(A + 27 ). Applying the quadratic formula, first to Eq. (A3), to obtain x then 
to A^ + 27 A — x = 0 gives the spectrum of . 

Ai = -7 ± y^72 + Xi (A4a) 




(A4b) 


A^^^ffA' 
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where ai = 7 ^ + (1 — p^) 7 i 72 + and oq = 7 i 72 (A^ + (1 — are the coefficients 

of the and terms of Eq. (A3) respectively. Simplifying Eq. (A4b) and substitnting 
Ap = -y/A^ + (1 — p 2 ) 7 i 72 yields Eq. (7) in the main text: 


Aj = -7 ± 


\ 


(7^ - AJ) ± ^(V-A?P + A 2 ( 7 i- 72 )! 


(A 6 ) 


Equation (A5) can be rewritten in the following convenient form 


Aj = -7 ± ApV’C^ - 1' 


1±^1 + T ]^ 


(A 6 ) 


where C = ^ cind 7 = are parameters defined in Eqs. ( 8 b) and ( 8 c) of the main 

text respectively. 


2. Underdamped Regime C = <C 1 

ZAp 

Consider now a V-system in the underdamped regime with 1. In such a system, the 
constraint on C, is equivalent to 

A2 = A2 + (l-p2)^i72»f (A7) 

In Eq. (A7), 7172 is the square of the geometric mean while 7 ^ is the square of the 
arithmetic mean. The geometric mean is less than or equal to the arithmetic mean and so 
7i72 a 7 ^- Noting that 0 < p < 1, Eq. (A7) gives the conditions for the underdamped 
regime 




(A 8 a) 


Ap ^ A 

(A 8 b) 

We now consider 7 = 

Applying the conditions in Eq. (A 8 b), we get 




(A9) 


We can now apply Eq. (A9) and the C ^ 1 condition to approximate Eq. (A 6 ) term by 
term using the binomial approximation. This gives: 
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(AlOa) 




2 


) 


i 


'i± VTT 


7 ]^ 



(AlOb) 


Substituting Eq. (AlOb) into Eq. (A6) gives the spectrum, {A*}, of A in the underdamped 
regime 


Ai = -7i 

(Alla) 

-^2 = “72 

(Allb) 

As, 4 = -7 ± 

(Allc) 

Now we proceed to determine the normal modes of (Eq. (Al)). 

E — \il), can be found with the usual method to give 

The eigenvectors. 

n n n 

v.cc [1,0,0, ] 

(Al2a) 

V2 0C [0,1,0, ] 

(Al2b) 

Vg OC [0,0, 1, 1] 

(Al2c) 

V4 OC [0,0,1, -1] 

(Al2d) 

We now use the general variation of parameters method [42] with the initial conditions, 

= 0, appropriate for excitation from the ground state to get the general form of the weak 


solutions of the master equations Eq. (lb) 
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(A13) 


pt 4 „t 

^ V / cis(vi • d) 
Jo io 


x(t) = e 


,Ai(t-s) 


Substituting Eqs. (Allc) and (Al2d) into Eq. (A13), evaluating the integrals gives, to 
lowest contributing order in ( the density matrix elements 


61 , 62 ' 


pil,e2it) = 


P^/rW2 


{t) = n{l — e 

(Al4a) 

Py/n^ -7t • /A 
= e ^ sm(Apt) 

L\p 

(Al4b) 

p '71^ _ p 72^7 

^'(cos(Apt) 1 ) 2 j 

(Al4c) 


Equation (A14c) gives the result presented in the main text. 

3. Overdamped Regime C 1 

Consider now the overdamped case, characterized by ^ 3 > 1. This is equivalent to 7 ^ Ap. 
As in the underdamped case we consider the parameter rj 


V = ^ 


- 72I _ A|7i - 72I ^ 2Ap 


1-2 A2| - ^2 <^«1 (A15) 

\T - ^p\ T 7 

We now apply a term by term binomial approximation to Eq. (A 6 ), this time using 






2(2 


(A16) 


Substituting the binomial approximations in ( ^ (Eq. (A16)) and in 7 (Eq. (AlOb)) into 
Eq. (A 6 ) we obtain the spectrum of in the overdamped region. 


Ai = -27 
A 2 


Ao — —- 


27 


(Al7a) 

(Al7b) 
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A2,3 = -7(1 ±*|) ~ -7 


(Al7c) 


Proceeding to find the eigenvectors of corresponding to the spectrum [Eq. (Al7c)] 
through the standard method: 


vi oc [rl,r2,p^/r^,0] 


(Al 8 a) 


V2 oc [r2,ri, -pv'r^,0] 


(Al 8 b) 


V3 oc [0,0,0,1] 


(Al 8 c) 



(Al 8 d) 


We note here that vi = d (Eq. (Al 8 a)) indicating that incoherent pumping drives V-type 
systems into the statistical mixture pd discussed in the main text (Eq. (3)). In contrast, V 2 
(Eq. (A18b)) represents the decay from pd to the rate-law predicted mixture peq (Eq. (5)). 

Finally, substituting Eqs. (Al7c) and (Al 8 d) into Eq. (A13) and computing the integrals 
we obtain the density matrix elements 



(Al9a) 



(Al9b) 


as discussed in the main text. 

4. Critically Damped Regime C = Ij 7 = = 0 

Consider the final regime discussed in the main text, the critical regime (C = 1 and 
7] = 0). In this case, the spectrum, {A*}, (Eq. (A 6 )) collapses to a single eigenvalue A = — 7 . 

Note that 7 = 0 can occur in two cases, if A = 0 or if 71 = 72 . Consider first, the case 
of A = 0 . In this case C = 1 (and so Ap = 7 ) implies that p = 1 and 71 = 72 . Substituting 
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this back into the equation Ap = 7 shows that this is simply the trivial case of uncoupled 
states (71 = 0 = 72). Consider next, the 71 = 7 = 72 case. In this case, C = 1 implies that 
A = p7. 

In the critically damped regime, the matrix becomes defective and cannot be diago¬ 
nalized. The variation of parameters procedure still provides the solution to the Quantum 
master Equations Eq. (lb) in the form. 

x(f) = [ (A20) 

Jo 

Where the exponential of A can be found by putting A into Jordan Canonical Form. This 
can be done by determining the generalized eigenvectors. Proceeding to find the generalized 
eigenvectors of rank j, using the standard procedure we obtain 






J2) 


J3) 


oc [ 1 ,- 1 , 0 , 0 ] 

(A 21 a) 

cx [ 1 , 1 , 0 ,- 1 ] 

(A 21 b) 

oc [ 0 , 0 , A-\ 0 ] 

(A 21 c) 

oc [ 0 , 0 , 0 , A- 2 ] 

(A 21 d) 


In the basis given by Eq. (A21d) the coefficient matrix takes the Jordan Canonical Form 


A = 


(A 22 ) 


/-7 0 0 0 \ 

0 -7 1 0 

0 0 -7 1 

y 0 0 0 —7 

Finally, substituting Eqs. (A22) and (A21d) into Eq. (A20) and evaluating the integrals 
we obtain the dynamics of a V-system in the critical regime 


pei,eA^) =n(l-e 


(A23a) 




(A23b) 
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(A23c) 


Piue2i^) = 


this reproduces the dynamics discussed in the main text. 


[1] R. E. Blankenship, Molecular Mechanisms of Photosynthesis (Blackwell Science, Oxford, 2002) 

[2] Y. C. Cheng and G. R. Fleming, Annu. Rev. Phys. Chem. 60, 241 (2009). 

[3] P. Wiirfel, Physics of Solar Cells: From Principles to New Concepts (Wiley-VCH, Weinheim, 
2005) 

[4] F. Rieke and D. A. Baylor, Rev. Mod. Phys. 70, 1027 (1998). 

[5] D. Segal and A. Nitzan, J. Chem. Phys. 122, 194704 (2005); D. Segal, Phys. Rev. B 73 , 205415 
(2006); 

[6] C. Creatore, M. A. Parker, S. Emmott, and A. W. Chin, Phys. Rev. Lett. Ill, 253601 (2013). 

[7] Y. Zhang, S. Oh, F. H. Alharbi, G. S. Engel, and S. Kais, Phys. Chem. Chem. Phys. 17 , 5743 
(2015). 

[8] M. O. Scully, K. R. Chapin, K. E. Dorfman, M. B. Kim, and A. Svidzinsky, Proc. Natl. Acad. 
Sci. USA 108, 15097 (2011). 

[9] K. E. Dorfman, D. V. Voronine, S. Mukamel, and M. O. Scully, Proc. Natl. Acad. Sci. USA 
110 , 2746 (2013). 

[10] G. D. Scholes, G. R. Fleming, A. Olaya-Castro, and R. van Grondelle, Nat. Chem. 3 , 763 

( 2011 ). 

[11] D. M. Jonas, Annu. Rev. Phys. Chem. 54, 425 (2003). 

[12] G. S. Engel, T. R. Calhoun, E. L. Read, T. K. Ahn, T. Mancal, Y. C. Cheng, R. E. Blankenship, 
and G. R. Fleming, Nature (London) 446, 782 (2007). 

[13] E. Collini, C. Y. Wong, K. E. Wilk, P. M. Curmi, P. Brumer, and G. D. Scholes, Nature 
(London) 463 , 644 (2010). 

[14] A. Chenu and G. D. Scholes, Annu. Rev. Phys. Chem. 66, 69 (2015). 

[15] N. Lambert, Y.-N. Chen, Y.-C. Cheng, C.-M. Li, G.-Y. Chen, and F. Nori, Nat. Phys. 9, 10 
(2013). 

[16] V. Tiwari, W. K. Peters, and D. M. Jonas, Proc. Natl. Acad. Sci. USA 110, 1203 (2013). 


26 




[17] A. Chenu, A. M. Branczyk, G. D. Scholes, and J. E. Sipe, Phys. Rev. Lett. 114, 213601 (2015). 

[18] L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, 
Cambridge, UK, 1995), Chap. 13. 

[19] X-P. Jiang and P. Brumer, J. Chem. Phys. 94, 5833 (1991). P. Brumer and M. Shapiro, Proc. 
Natl. Acad. Sci. USA 109, 19575 (2012). 

[20] 1. Kassal, J. Yuen-Zhou, and S. Rahimi-Keshari, J. Phys. Chem. Lett. 4, 362 (2013). 

[21] T. Grinev and P. Brumer, J. Chem. Phys. 143, 244313 (2015). 

[22] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, Atom - Photon Interactions: Basic 
Process and Applications (Wiley-VCH, Weinheim, 2004). 

[23] R. Loudon, The Quantum Theory of Light, 3rd ed. (Oxford University Press, Oxford, 2000). 

[24] M. Fleischauer, C. H. Keitel, M. O. Scully, and C. Su, Opt. Commun. 87, 109 (1992). 

[25] V. V. Kozlov, Y. Rostovtsev, and M. O. Scully, Phys. Rev. A 74, 063829 (2006). 

[26] D. Gelbwaser-Klimovsky, W. Niedenzu, P. Brumer and G. Kurizki, Sci. Rept. (in press) 

[27] M. Kiffner, M. Macovei, J. Evers, and C. H. Keitel, Prog. Opt. 55 , 85 (2010), sec. 3.2.5. 

[28] G. S. Agarwal, Quantum Statistical Theories of Spontaneous Emission and their Relation to 
Other Approaches (Springer-Verlag Berlin, 1974). 

[29] M. Fleischauer, A. Imamoglu, and J. P. Marangos, Rev. Mod. Phys. 77, 633 (2005). 

[30] S. Harris, Phys. Today 50 , 36 (1997). 

[31] L. V. Hau, S. Harris, Z. Dutton, and C. H. Behroozi, Nature 397 , 594 (1999). 

[32] M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge University Press, Cambridge, 
UK, 1997). 

[33] A. Dodin, T. V. Tscherbul and P. Brumer, (manuscript in preparation). 

[34] Note that while Fig. 3 of Ref. [36] shows the absolute magnitude of the coherence vs. time, the 
time evolution is not discussed in that paper, which focuses on the coherences in the long-time 
limit. 

[35] B.-Q. Ou, L.-M. Liang, and C.-Z. Li, Opt. Commun. 281, 4940 (2008). 

[36] G. S. Agarwal and S. Menon, Phys. Rev. A 63, 023818 (2001). 

[37] G. C. Hegerfeldt and M. B. Plenio, Phys. Rev. A 47, 2186 (1993). 

[38] T. P. Altenmiiller, Z. Phys. D 34, 157 (1995). 

[39] S. Menon and G. S. Agarwal, arXiv:quant-ph/9902021vl (1999). 

[40] T. V. Tscherbul and P. Brumer, Phys. Rev. Lett. 113 , 113601 (2014). 

[41] T. V. Tscherbul and P. Brumer, J. Chem. Phys. 142, 104107 (2015). 


27 


[42] W. E. Boyce and R. C. DiPrima, Elementary Differential Equations, 9-th Edition (Wiley, NY, 
2008) 

[43] D. Morin, Introduction to Classical Mechanics (Cambridge University Press, UK, 2008). 

[44] Z. S. Sadeq and P. Brumer, J. Chem. Phys. 140, 074104 (2014). 

[45] T. V. Tscherbul and P. Brumer, Phys. Rev. A 89, 013423 (2014). 

[46] J. Olsina, A. G. Dijkstra, C. Wang, and J. Cao, arXiv:1408.5385vl (2014). 


28 


